Landau-Zener transition of a two-level system driven by spin chains near their critical 

points 
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The Landau-Zener(LZ) transition of a two-level system coupling to spin chains near their critical 
points is studied in this paper. Two kinds of spin chains, the Ising spin chain and XY spin chain, 
are considered. We calculate and analyze the effects of system-chain coupling on the LZ transition. 
A relation between the LZ transition and the critical points of the spin chain is established. These 
results suggest that LZ transitions may serve as the witnesses of criticality of the spin chain. This 
may provide a new way to study quantum phase transitions as well as LZ transitions. 
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INTRODUCTION 

Quantum information processing promotes renewed at- 
tentions in quantum two-level systems in recent years. A 
number of two-level systems have been tested as good 
candidates of qubits that are the least units in quan- 
tum information processing A good qubit requires 
that the qubit is well isolated from its environments and 
easy to manipulate. There are many ways to manipulate 
qubits, one of them is to use Landau-Zener (LZ) sweeps, 
which have been realized in recent experiments on super- 
conducting qubits @, H, HJ] . 

The LZ transition occurs when two instantaneous 
eigenvalues of a quantum system come close together due 
to the parameter change. It has attracted attentions for 
decades since the work in the early 1930s in slow atomic 
collisions [a, |g, |7fl and spin dynamics [8( . The LZ theory 
has found many applications, such as the mentioned ma- 
nipulation of qubits, the enhancement of the macroscopic 
quantum tunneling |9j and the control of the transition 
probability as well as the quantum phase factor [n| ■ Re- 
cently, the LZ effect was proved to be useful in quantum 
information processing such as the preparation of quan- 
tum states in circuit QED systems [111 ], creating entan- 
gled modes in cavities [12| and fault-tolerant single-qubit 
gate operations [l3T | . In practice, the qubit is always influ- 
enced by its environment, leading to decoherence in that 
system. Therefore, taking the environment into account 
when study the LZ transition is a practical considera- 
tion. The LZ transitions in two-state systems dissipa- 
tively coupled to their environments have been studied 
extensively 



14, 15, la 17, hj, those results show that 



the two-level system is robust against the influence of 
environment. Those studies do not consider the interac- 
tions among the particles in the environment, hence the 
information of the particle-particle correlation can not 
reflected in the LZ transitions. It is well known that the 
spin-spin coupling in the spin-chain may result in quan- 
tum phase transitions, due to the quantum fluctuation at 
zero temperature (l9| . One dimension spin chains which 



may be solved exactly are rather attractive in such stud- 
ies, especially by using the concepts developed in quan- 
tum information theory, such as entanglement of spins in 
spin chains 12011 , as well as the geometric phase 21 1 and 
the fidelity [22| . Other relations have already been estab- 
lished between the quantum phase transitions in environ- 
ment and the decoherence of system [23j, entanglement 
24 1 , geometric phase 2iJ 26| , where the property of sys- 
tem play the role of detector of the environment. These 
studies have also stimulated us to consider LZ transi- 
tions in the environment of spin chains, and such a study 
may provide a new method to investigate quantum phase 
transitions. Actually, the relation between LZ transitions 
and the quantum phase transitions has been considered 
in Ref . [27j , where they studied the LZ transitions of the 
spin chain itself under a time-dependent evolution, here 
we shall consider the LZ transition of a coupled time- 
dependent two-state system and regard the chain as an 
environment. 

In this paper, we will investigate the LZ transitions in 
two-state systems surrounding by spin chains in trans- 
verse fields near their critical points. We will focus on 
two kinds of spin-chain, one is the Ising spin chain and 
another is XY spin chain. As we shall show you, the envi- 
ronment's properties can also be reflected in the systems' 
LZ transitions. 



ISING SPIN CHAIN AS AN ENVIRONMENT 

Taking a spin chain described by the Ising model as 
the environment, we restrict ourself to consider the case 
where the chain-system coupling only affects the spin flip 
of the two level system. The Hamiltonian that governs 
the dynamics of the whole system (two-level system and 
the chain) reads, 
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Here 



is the standard Landau-Zencr Hamil- 



tonian! with a z = |T> (Tl - II) (41, and a x = |t)(|| + ||)(t|, 
|t), 1 1) stand for the excited and ground state of the two 
level system, and this part of the Hamiltonian describes 
the situation where the dynamics is restricted to two lev- 
els that are coupled by a constant tunnel matrix clement 
and cross at a constant velocity, a x , and cr| stand 
for the Pauli operators of the spin at site j in the Ising 
chain; J and A describe the strength of the spin coupling 
and the transverse field. The periodic condition has been 
chosen to the spin chain, and we choose M = (N — l)/2 
for odd N, where N is the number of the spins in the 
chain, g is the coupling coefficient between the two-level 
system and the surrounding Ising chain. 

It is convenience to continue the calculation in the in- 
teraction picture, in order to calculate the probability of 
the qubit state flips due to the LZ sweep, and we may di- 
vide the Hamiltonian into two parts as H(t) = Ho(t) + V, 
where H (t) = f a z - JJ2 3 ( cr j (7 j+i + Xa V is the bit - fli P 
free Hamiltonian and V = cr x (~ — dJ2j a j) describe the 
bit-flip interaction. 

Obviously, H (t) may be diagonalized as long as we 
diagonalize the part of Ising spin chain Hamiltonian, and 
this may be realized through Jordan- Wigner transfor- 
mation ai = Yl m <i a m(~ a i + i <J i)/2i Fourier transfor- 
mation, Cfc = ai exp(—i2nlk/N), and Bogoliubov 



transformation, bk = Ck cos -# 



J 



sin t£ , and in the 



diagonalize procedure, we have defined cos Ok = £/c/£fcj in 
which £k and £fc arc defined as £k = 2J(A — cos ^jp-), while 



N 

£ fc are defined as £ k = 2 J [cos ^ - A] 2 + sin 2 
These transformations have transformed the spin oper- 
ators into the quasi fermion operators in the momentum 
space, which may greatly simplify our following studies. 
From this procedure, we finally obtain the diagonalized 
expression, Ho(t) = ^^ z +^2k^bkbk, following which we 

will get the expression of Uo(t) = exp[— i f_ H (t')dt'}, 
i.e., 



U (t) = exp(-i Y] ^-b\b k t) ■ exp(--^-cr 
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The Hamiltonian in the interaction picture may be ob- 
tained through H(t) = U (tfVU (t), and we may rewrite 
V into the form described by b\ and bk, which has been 
defined^ in the above diagonalize procedure, in order to 
bring H(t) into a useful form for our following study. De- 
fine J x — a x , which is a magnetic moment operator 
of the chain, and V = <j x W in which W is defined as 
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Define the product state |n) = \ti-m) ■■■\no)\rii) ...\um) , 
which stands for the state of the quasi fermion envi- 
ronments (Ising chain), n = n , n 1; n 2 , %} 
with m = 0, 1, f2 = {uj-m, uj , u>i, uj 2 , um} with 
^k = £k/h> and we may change the sum J2k n k£,k/'h into 
the inner product form n • O for simplicity. 

Then, using @, © and the completeness of I = 
J2 n l n )( n L we m ay obtain H(t) easily, 

H(t) = ^exp[i(m - n) • ftt}W mn \m) (n\ 



,ivt 



ivt 



[exp( — )|T)(||+exp(- — )U)(T|], (4) 



in which we have defined W mn — (ml Win) 



LANDAU-ZENER TRANSITION PROBABILITY 

Now we can calculate the LZ transition probability 
based on the above preparing work. Suppose that at time 
t = —oo the two-level system is in its excited state |f) and 
the Ising spin chain system starts in its ground state, i.e., 

|0) = rL>o( cos \ + * sin ^4 c -fc)l )fcl°)-fc which satis- 
fies b±fc|0) = 0, and the state of the whole system may be 
expressed as |V>(— oo)) = ||)|0). We now begin to calcu- 
late the survival probability of the initial state f) at time 
t = oo, i.e., P|^|(oo)=|(t|'0(oo))| 2 . The evolution oper- 
ator t/(oo,-oc) with U(t 2 ,ti) = T exp[-i drH(r)], 
may be expressed into a time-ordered expansion, with 
ti < t% < ... < <2fc-i < t 2k in the interval (— oo, oo), and 
only the even powers of H(t) will contribute to P|— >f (oo). 
The perturbation series for (f \ip (oo)} can be expressed as, 



v -| />oo />c 

^ Vt / n (2fc) n (2fc-l)I / l /r n (2fe-l) n (2fc-2) ...W n (l) 



dt 2 — I dt 2 k 
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in which we have defined w\ — n' 1 ^ • f2, wi — (n^ — 
n'' -1 )) ■ O, (2 < I < 2k) for simplicity, and n' 1 ' denote 
the state of the chain after the i-th interactions. 

In order to work out the above integration, it is ad- 
vantageous to make the variable transformations. Intro- 
duce a set of new variables as [ljjj, x q = Y^a=~i 1 (~ ^) l+1 ti, 



Vq = 

El-i 
q=l Vq 



T-2q 



^2(7-1, and consequently xi 



t 2 i- 



Then the above perturbation series 



(3) ([5]) for (1"|V>(oo)) can be changed to the form Eq. 
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An analogous analysis to references EI El will help us 
to work out the above integration. Under permutation of 
the xi, the integrand is not symmetric. If we transform 
the x\ into new variables s\ = x\ and si = xi — Xi—i 
for I = 2,3, ...,k, then xi — J2q'=i s q'- ^ is easy to 
find that the dsi-integral will give us the (5-function 
{27T%/v)6(Y J i =1 Vi + |(^2/ + w 2 i-i) )■ And since the ini- 
tial state is |T)|0), we can obtain Y^i=i( w 2l + w 2i-i) = 
n' 2fe )-f2 > 0, further more, the variables yi > 0, therefore, 
the integral on delta- function can only exist in the sub- 
space yi=y2=- ■ ■=J/fe=0 and only in the condition that 
the vector n( 2fe )=0. Then we can return to the above in- 
tegration with the condition of yi = 0, (I = 1, k), and 
since within this subspace the integrand symmetric in the 
variables Xi, we may symmetrize the a;;-integrals, by re- 
placing the integrals of xi to (1/A;!) /_ dx\ . . . J_ dXk- 
After the Xi integrals, then the ^/-integrals can also be 
evaluated as well by using the property of (5-functions, 
J 5(y)dy = 1/2. From the above time integrals we have 
noticed that when the environment starts in the ground 
state |t) |0), only the (2fc)th-order processes contribute 
to the survival probability Pj_^(oo), and they should 
satisfy n( 2 )=n( 4 )= . . . =n( 2fc )=0, then the quasi fermion 
environments will end up in their initial state |0) in case 
the two-level system ends up in |f). However, the time 
integrals of the equation do not prohibit the occupation 
of the states |T)|n ^ 0) at intermediate times, and they 
also do not restrict the intermediate environment states 
| n (2/-i)^ the vanishing matrix elements W ma 

give the further restrictions to the integral. 

According to the above analysis and calcula- 
tions, we find that (|| - i/'(oo)) can be simplified into 
cxp(— 7rr 2 /?w)|0}, in which the parameter T 2 have been 
defined as T 2 = *£ n W 0ll W n o = (0\W 2 \0) and we may 
obtain the exactly LZ transition probability for a qubit 
coupled to the Ising spin-chain in a transverse field, 



-2tiT 2 /hv 



P T ^;(oo) = l-P T ^ T (oo) = l-e 



in which the parameter T 2 can be work out as 



(7) 



4 -g(J x )of + g 2 (A,r i 
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= ( T -sE cos ^) 2 +3 2 E sin2 ^' (8) 
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FIG. 1: (color online) {J x )o and (AJ X ) 2 (solid line) as well as 
their derivatives by A (dot-star line) as the functions of the 
parameter A. We set N = 200 in the numerical calculation. 



where the term {J x )o = J2k>o cos @ k * s J us ^ ^ ne ex P ec t- 
ing values of the magnetic moment of the spin chain, 

and AJ X ee V((^) 2 )o - {J x )l = ^Ek>o s ' m2 ° k is thc 
variance of the magnetic moment at the ground state 
0). In case of the coupling coefficient 5 = 0, the ex- 
pression is in consist with the usual result Pj-_>|(oo) = 
1 — exp(— irA 2 /2Tiv), and when A = 0, the transition are 
completely determined by thc environment chain. 

The final expressions (J7J) and ([5]) tell us that the LZ 
transition of the central two-level system depends on 
both the expecting value of the magnetic moment J x and 
its variance at ground state of the spin chain. Both the 
expecting value and its variance are closely related to the 
strength of the transverse field, which will affect the LZ 
transition equation through 9k, thus we may conjecture 
that the message of quantum phase transition happens 
at the certain strength of the transverse field may also be 
reflected on the LZ transition probability of the two-level 
system. 

Fig. [T] shows (J x )o and (AJ X ) 2 as functions of the 
transverse field strength A, as well as their derivatives 
with respect to A. When the strength of the transverse 
field is in the region 0<A<1, only the expectation value of 
J x changes, while the variance part of T 2 does not. This 
feature tells us that, the change of LZ transitions mainly 
caused by the expectation value of magnetic moment in 
this case; however, when A>1, the variance become de- 
creased with the transverse field, both of them will affect 
the LZ transition and the expectation values of magnetic 
moment may become the mainly causation. The expec- 
tation value changes sharply when A — * 1, indicating that 
their derivatives may reveal the singularity near the crit- 
ical points perfectly. The above property becomes more 



FIG. 2: (color online) In units of J 2 , F 2 and its derivative 
<9F 2 /<9A as the functions of A and A . We set N = 200, 
g = 0.1J, A m ax = 20 J, in the numerical calculation. 
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FIG. 4: (color online) In the environment of XY spin chain, 
LZ transition probability P and its derivative dP/dX as the 
function of A and 7. We set N = 200, g = 0.1J, A = 5J, 
v = 50J 2 /h, in the numerical calculation. The panel shows 
dP/dX changes with A in case of 7 = 0.2, 0.5, 0.8, 1.0. 



FIG. 3: (color online) LZ transition probability P and its 
derivative dP/dX as the function of A and A. We set N = 

200, g = 0.1J, A max = 20J, v = 50J 2 /h in the numerical THE XY SPIN CHAIN AS THE ENVIRONMENT 

calculation. 



clear when N increases, this result comes from the analy- 
sis given m In the following we shall present 
the results with a specific number N = 200. 

Since the magnetic moment of the spin chain changes 
fast at the critical points of the chain, then we may ex- 
pected that the LZ transition can also reflect such prop- 
erty of the spin chain by the expression of (O and |8]). 
Fig. [2] shows the function of T 2 as well as its derivative 
relation with A and A, we set the number of the chain 
N = 200, g = 0.1J in the numerical illustration. It is evi- 
dent that when the strength of the transverse field A — * 1, 
r 2 also changes sharply. As mentioned above, in order 
to reveal this effect of quantum phase transitions more 
clearly, we may pay more attention to the derivative of 
r 2 by A, which perfectly shows the quantum critical phe- 
nomenon. Just as we have expected, this property has 
been inherited very well by the LZ transition probabil- 
ity Pf^x(oo) and its derivative, as Fig. [3] shows, where 
we have set v = 50 J 2 /ft in the calculation. It is in evi- 
dence that the critical point is reflected perfectly well in 
the derivative of LZ transition, which is consist with the 
above analysis. 

Different values of A also affect the property of LZ 
transition probability. From Eq. © and Fig. [3J we 
may find that when A > 0, the transition probability de- 
clines first, and then reveals again, however, the sharply 
changed location is not influenced, which can also give us 
a good reflection of quantum phase transitions. When A 
is large enough, the transition probability will not reveal 
again in the case of A > 1. 



Now we consider the case where the XY spin chain acts 
as the environment. Straightforward calculation shows 
that to get the result in this situation, we need to re- 
place -JEjLjfi^^+l + ^crjo-j+i + Act- ) with 

~JJ2jL-M( a j cr j+i + AcrJ) in E< F ©■ Here 7 measures 
the anisotropy in XY spin-chain. The XY Hamiltonian 
will turn into the transverse Ising chain for 7 = 1, and 
the XX chain in transverse field for 7 = 0. 

By the same procedure, we can obtain the LZ tran- 
sition probability in this case, the equations are noth- 
ing but changing the definition of 8 k and £/. in the 
above discussions, i.e., cos 9k = £fc/£fc, in which Ek = 



2J(A 



2?rfc ' 
N 



and £k are now defined as 



2.7 



cos 



2-rrfc 
JV 



A] 2 



7 2 sin 2M 



Following the same calculation, we then obtain the 
same expression as and ([5]). Fig. 0] shows the relation 
between LZ transition probability P and the anisotropy 
parameter, as well as its derivative dP/dX. We set 
N = 200, g = 0.1J, A = 5J, v = 50J 2 /^, in the nu- 
merical calculation. It shows that when 7 changes from 
1 to 0, the critical behavior is reflected quite well in the 
LZ transition, and the critical line A = 1 is clearly re- 
flected in the derivative of the LZ transition probability. 
We may deduce that the distinctive feature of the XX 
chain environment is a sharp discontinuity in the deriva- 
tive dP/dX at A = 1, this is reflected more clearly from 
the panel in Fig. |U This result confirm our above predic- 
tion that the LZ transition can reflect the critical points 
of the environment. 
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CONCLUSION 

In summary, we have studied the Landau-Zener transi- 
tions in two-state systems coupling to the Ising spin chain 
and XY spin chain in transverse fields. We have calcu- 
lated the exact expressions of the LZ transition proba- 
bilities of the two-state systems and analyzed the rela- 
tion between their properties and the occurrence of the 
quantum phase transitions in the chain. The results show 
that the LZ transition are determined by the spin chains' 
magnetic moments and their variance. As the magnetic 
moments of the chains contain the information of quan- 
tum phase transitions, the LZ transitions may act as the 
witnesses of quantum phase transitions in the chains. 
Our results suggest a rather intriguing relationship be- 
tween LZ transitions and the environments' properties, 
and therefore the results may provide a new way to study 
the phenomenon of quantum phase transition as well as 
Landau-Zener transition. 

This work was supported by NCET of M.O.E, and NSF 
under grant No. 60578014. 
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